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SUMMARY 

The thrust, houndary-layer, and heat-transfer characteristics were 
computed for nozzles having radial flow in the divergent part. The work- 
ing medium was air in chemical equilibrium, and the boundary layer was 
assumed to be all turbulent. Stagnation pressure was varied from 1 to 
32 atmospheres, stagnation temperature from 1000° to 6000° R, and wall 
temperature from 1000° to 3000° R. Design pressure ratio was varied 
from 5 to 320, and operating pressure ratio was varied from 0.25 to 8 
times the design pressure ratio. 

Results were generalized independent of divergence angle and were 
also generalized independent of stagnation pressure in the temperature 
range of 1000° to 3000° R. A means of determining the aerodynamic ally 
optimum wall angle is provided. 


INTRODUCTION 

In general, the most important characteristic of an exhaust nozzle 
for a propulsion device is the amount of thrust that it supplies at a 
given pressure ratio, stagnation temperature, and stagnation pressure. 

This number is usually compared with the performance of an ideal nozzle; 
that is, a nozzle that isentropically expands the exhaust gas to free- 
stream pressure and a purely axial direction at its exit station. Cal- 
culations to determine the ratio of actual to ideal thrust generally de- 
pend on pressure or temperature level of the process even if viscosity 
is neglected; this is due to changes in specific heat and molecular 
weight of the working fluid. At moderately high temperatures the amount 
of heat transfer becomes important for structural considerations so that 
the boundary layer, which in turn determines the heat-transfer coefficient, 
can no longer be neglected. Since nozzle exit angle and design pressure 
ratio (or area ratio) also must be considered, the selection of an opti- 
mum nozzle for a particular application becomes a sizeable task. 
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The figures of this report were prepared to reduce repetitious com- 
putation of thrust and heat-transfer characteristics of radial-flow noz- 
zles using air as a working medium. The general trends and order of 
magnitude of viscous effects may aid in the selection of nozzles using 
fluids other than air. Boundary-layer momentum thickness and form factor 
are also included as a possible aid to those engaged in developing sepa- 
ration triggering devices. Radial-flow nozzles are defined herein as 
nozzles that have all the potential-flow velocity vectors of the diver- 
gent section directed radially outward as if originating from a point 
source located ahead of the nozzle throat. All the potential-flow prop- 
erties are therefore constant on any spherical surface with center at 
the hypothetical point source. The bounding walls of the divergent sec- 
tions are conical. The convergent sections considered are those plotted 
in figure 1. 

All calculations are based on air at chemical equilibrium with a 
composition of 20.95 percent oxygen, 78.09 percent nitrogen, 0.93 per- 
cent argon, and 0.03 percent carbon dioxide at standard conditions. Cal- 
culations were made for nozzles having design pressure ratios of 5 to 
320 and inlet stagnation temperatures and pressures from 1000° to 6000° R 
and 1 to 32 atmospheres, respectively. Although inlet stagnation temper- 
ature of 1000° R is well below that of much practical application, it is 
in the range of readily available laboratory air supply and therefore in 
the range of considerable experimentation. Comparison of the low- and 
high-temperature calculation may then aid in the utilization of low- 
temperature data to high-temperature applications. Nozzle wall temper- 
ature was held constant at 1000°, 2000°, and 3000° R. These constant 
wall temperatures correspond either to the case of a heavy-walled nozzle 
when the longitudinal transfer of heat can be very large compared with 
that flowing through the boundary layer or to the case of a nozzle having 
coolant flow distributed to hold wall temperature constant. The boundary 
layers calculated were assumed to be all turbulent. 


SmBOLS 

A area on a spherical surface (fig. l) 

a speed of sound 

B coefficient in eq. (3) 

F thrust in a direction parallel to noz 2 le axis 

g gravitational constant 

H profile parameter, 8 * / 9 
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enthalpy 

length 

Mach number 

rate of mass flow 

Prandtl number 

pressure 

rate of total heat transfer 

rate of heat transfer per unit area 

radial distance from hypothetical point source in radial flow 

distance normal to nozzle axis 

temperature 

velocity 

rate of weight flow 

transformed distance along nozzle wall 
distance along nozzle wall 
transformed distance normal to nozzle wall 
distance normal to nozzle wall 
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|jL absolute viscosity 

v kinematic viscosity 

a gas constant 

p mass density 

Subscripts : 
av average 

aw adiabatic wall 

b base 

e inside nozzle at exit station 

g expression independent of nozzle size and wall angles 

I ideal 

i incon^ressible or transformed 

opt optimum 

p potential stream 

ref reference 

s stagnation 

std standard conditions (for pressure, 1 E.tm) 

T total 

t potential-stream stagnation 

v viscous 

w wall 

00 undisturbed free-stream conditions exi emal of nozzle 

1 1 foot 
Superscript: 

* refers to nozzle throat except when used with displacement 

thickness 
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ANALYSIS 

The flow through the nozzles was considered in two parts: a poten- 

tial flow in the major part of the nozzle , and a boundary layer in which 
viscosity and conductivity were considered. The charts of reference 1 
were used to determine pressures, temperatures, and enthalpies of the 
air in the potential-flow region. The flow direction and static pres- 
sure normal to the flow direction were considered constant in the bound- 
ary layer, so that nozzle wall pressures were determined by the potential 
flow. 


The calculations for the development of the viscous layer of air 
along the walls were based on air with a constant molecular weight and 
essentially constant specific heat but with the proper boundary condi- 
tions at the edge of the potential layer > Some inherent inaccuracy in 
the upper temperature levels considered therefore may be present in the 
heat- transfer and boundary-layer data presented. 


Basic Equations 


Thrust ratio and boundary layer * - Ideal thrust coefficient Fj is 

defined as the thrust obtained by a nozzle expanding in all-isentropic 
flow to ambient pressure (p^) with all the terminating velocity vectors 
in an axial direction at the nozzle exit. Also, by definition, the actual 
and ideal nozzles must pass the same amount of air per unit of time. 

The total mass flow through the actual and/or ideal nozzle is 


m == 



dA 



( Pp U p - pU)dA 


where A is the area on a spherical surface with the center at the hy- 
pothetical point source of the radial flow and the subscript on A de- 
notes the boundary of A. 

Utilizing the definition of displacement thickness gives: 


m 



dA. - 2jtr v 6*p p Up 



dA 


where second-order terms of the boundary- layer parameters are neglected, 
so that 


r 6 
w 


(r w - &*)&* 
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The quantity Ac* is herein the flow area of a nozzle with all- 

potential flow, whereas the actual nozzle radius must he increased by 5 
to keep the rate of mass flow the same in both cases. 


The actual jet thrust can be written in integral form as 


-I 


K 


7ip = j pU cos a dA + 

0 -o 


rK 

J (p - pJ 


cos a dA 


where all terms except p^ are evaluated inside the nozzles at the exit 
station. It is convenient to separate Frp into potential- flow and 
viscous-flow components as follows: 


f t = 


g f *' w 2 

J P p^p COS a cLA + J p p% cos a ^ 

A & * 


0 


/ Aw 2 p 

(ppUp - pUUp + pUUp - pU )cos a dA + 

a 5 

/•Ag* /*h w 

/ (p - P<Jcos a dA + J (p - p^) cos a 


dA 


so that 


[mUp + (p - pjA & *]cos av a - 2jtr w [p p uj;e - (p - pJ&*]cos 


mU T 


= _ £V 


(1) 


The potential- flow calculations were always made normal to a spher- 
ical control surface and were resolved into the axial-flow direction by 
multiplying by the average cosine of the aiigle between the flow direction 
and the nozzle axis. The development in appendix A shows that 


cos 


a 


av 


1 + COS (Xy- 
2 


( 2 ) 


for a radial-flow nozzle having circular cross sections normal to the 
axis . 
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Momentum thickness was computed by combining equations (20) and (23) 
of reference 2 to give 

0.823 


6 = 


TL 




r* 

m \0. 732 

± P ) 1.2155 


0.01173 j 

W r “ 

(H 

A v t 

0.2155 J 

A .£+0.2155 1.2155 1 

:) p r * l 

D 

© J 


where 


B = 4.2 + 1 


.563^ - l) 


and T re £ is evaluated by ^- re f and reference 1, where 


(3) 


ref l/ n w 


hp = % + 1 j + °- 22 ' 


L aw 


fp 


- 1 


and 


h 


aw _ 1 + x - 1 




“p 2 

Form factor H =6/0 was estimated by the expression 


_ h t + hjfgv 

"p 

which is developed in appendix B. 

Substituting the continuity equation 

k *\2. 


m 


= - 5TP^ Up 


cos av a 


and equation (4) in equation (l) results in 

A c 


f t U p + g(p - Peo ) 
pv = m 


cos av a 


Atf h t + 1 . 280 ^ 

Up - g(p - p J -^-1 II 


U-r 




(4) 


20 

( _ 6 ») GOS av a cos °w ( 5 ) 
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Heat transfer . - Combining equations (27) and (28) of reference 2 
results in the following expression for rate of heat transfer per unit 
area: 


\0.732 


0.123e 


2 VT^J ( Pr )ref 5 PpV h aw “ *V) 

q = /n, \0. 268/ \0. 268 

ft) (^) 

The total rate of heat transfer from the nozzle is 


( 6 ) 


Q = 2n 



r w q 3x 


(?) 


Generalization of Results Independent of 
Nozzle Size and Wall Angle 


Thrust and boundary layer. - The expression 


i 5, _ u p + gtp - tt . 

Gos av a F I U I \ F I 

is independent of nozzle size (scale) and wall angle. 

It is desirable to generalize F y /F j also. Let and be 

the distances of boundary-layer travel along the convergent and divergent 
section, respectively, of the nozzle. Assume that the boundary layer is 
all turbulent and that it has zero thickness at the entrance to the con- 
vergent section. Dividing all terms in equation (3) by L = + L 2 re- 

sults in 




( 9 ) 
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From figure l(b), 


L 2 = 


L w 


L W 



sm 


% 


sm 


°V 


if radial flow is assumed at the throat. 


Consider the family of convergent nozzles having lengths 
Li oc r*/sin such as sketched in figure 1(b). The convergent section 
is stretched as Oy decreases , so that any increment of length AL be- 
fore stretching results in a length AL 1 = (sin a/sin a 1 )AL after 


stretching. Then, 
% = 15° and r w 
point, 


for a fixed A^A*, L oc 
* 


L w 


or arbitrarily using 


sm oty' 

= a 1-foot nozzle throat radius as a reference 


L oc 


r* sin 15° 
r* sin 


( 10 ) 


From equation (9), 


r^ sin \ 0,823 
r* sin 15 ° ) 


(H) 


is independent of nozzle size (scale) and wall angle. Use of expression 
(11) in equation (5) then shows that 



is independent of nozzle size and wall angle for the family of nozzles of 
figure l(b). 


Vail angle for optimum thrust . - For radial-flow nozzles, an increase 
of wall angle decreases the axial components of thrust but also decreases 
the thrust loss due to viscosity. An optimum wall angle for maximum total 
thrust F t is therefore obtainable. Combining expressions (2), (8), and 
(12) gives 
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fr *2 fv 

F I " F I + F I 




1 + COS 0^ 


*\0.177 


r l \ (sin 15°\ 

,r * \ sin % ) 


0.823 


(cos a v )(l + cos cXy ) 


'gV-v 


Differentiating Frp/Fj with respect to and setting the result equal 

to zero yields 


/ *\° 
fe 


- 1 ” fel 


sin' 


2.823 




^ (sin ISO) 0 ’ 823 [(sin2a tf>opt ) (0.177 + 1>177 cos a w>opt ) + 0.833(l + cos <V, 0pt )] 


(13) 


V 


MV 

Equation (13) is plotted in figure 2 so that, if I— ^ 1 

\ r w/ 

* 

known, the optimum wall angle can be determined. 

Heat transfer. - Use of expression (ll. in ( 6 ) indicates 


°- 177 \ii 


& 

l Ft 


■s is 


/ # . \ 0.221 

fl sin “v \ 


\r w sin 15° 


(14) 


is independent of nozzle size and wall angle. 

Using expressions (l0) and (14) in (7) gives generalized Q as 





x 


w / ( Tj sm ay 
i yf sin 150 


^ 0.221 I s j n 2.50 r v 

V s in °v rp 


( 15 ) 




11 


Generalization of Results Independent of Stagnation Pressure 

When stagnation temperature is in the range 1000° = Trp = 3000° R, 

pressure has no effect on specific heat. In this range, however, the 
ratio of specific heats decreases with increase of temperature because of 
increases in the amount of the total internal energy content of the air 
molecules going into vibrational energy. Molecular weight is constant 
throughout this range. When temperature is increased much beyond 3200° R, 
molecular weight and specific heat both vary with pressure as well as with 
temperature (ref. l). 


Thrust and boundary layer . - Potential -flow thrust ratio can be 
written in terms of pressure ratio and thus is independent of stagnation 
pressure p^ for 1000° = = 3000° R. With equation (l) , 


F p U p + g(p - pj - 

f t u t 


cos av a 


or 


£e i 


F I cos av a 


= 1 E = 


I£ m! + (z- - 5=) 

Pt ^ \pt Pt / 




_ 

Pt , 

U T 

* w -n t\ 


In this lower temperature range, pressure enters into the momentum- 
thickness equation (eq. (3)) only through the kinematic viscosity term 
and can be factored out as: 


vv 


0. 01173 (oT t ) 


0.2155 


0.2155 


m -c 


Therefore, 


2155 1.2155 
r v 


e oc 



VO. 732 




r B \ T ref 


fit] 

\ T P/ 


1.2155 . 
3.268— r w 


0.823 


fc)' 


0.177 


(k) 


0.177 


( 16 ) 


and combining with expression (ll) results in 


0 (Pt.) 


0.177 (*1 sin a w 


\ r w sin 150 j 


0.823 
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being independent of wall angle, sizes, and pressure in the lower temper- 
ature range. 


Profile parameter is independent of pressure in this temperature 
range. Pressure enters into the viscous-thrust-loss expression, there- 
fore, only as 


' ± vO.177 

w 


since 


[ppUp - (p - pJ h ] 20 

0 


5 *Vp U I 


ccs av a cos % 




20 


(r, 


w 


- 5*) 


Pt 


MpMp 


J25 

\r t 

1 p p 


cos av a cos “v 


and 



Using equation (12) results in the expression 


/ \ . „ ,0.177 0.823 

( F v\ 0.177 F v ^ r vPt) sin “v 

\F t ) p t " F t / . 0 x 0. 822 / W x 

\ 1/ I (sin 15 u ) (cos av a)(.cos o^) 

o 

which is independent of nozzle size and wall angles as well as pressure 
in the lower temperature range (1000° to 3000° R). 

Heat transfer . - Using expression (16) and the equation of state 
in equations ( 6) and (7) shows that 

0.779 
q °c p t 


and 


0.779 
Qoc p t 


in the lower temperature range. 
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Therefore, from expressions (14) and (15), 


' * \ 0.221 


-0.779 

Pt 


and 

( . \0. 779 

^ Sln °v' -0.779 

4 *1.779 p t 

r 
w 

are Independent of nozzle size, wall angle, and stagnation pressure in 
the lower temperature range. In the presentation of the lower temper- 
ature range computations, p^ is arbitrarily divided by 2p s td raa ^ e 

it dimensionless in these parameters. 


DISCUSSION OF RESULTS 


Thrust 


Potential flow . - The thrust ratio obtainable from potential flow 
through radial-flow nozzles was computed for nozzle design pressure ra- 
tios of 5, 10, 20, 40, 80, 160, and 320. Thrust ratio is obtained from 
figure 3 by multiplying the ordinate by the cos av a of the nozzle being 

considered. This average cosine of the flow direction can be obtained 
from equation (2). 


The thrust curves are characterized by their rapid drop as the oper- 
ating pressure ratio is decreased below design pressure ratio, 

p+ £t 

Poo 


_t 

Poo 


© 


Pt 

p* 


= — < 1, and by their relatively constant value as 

Pro 


design 


P e /Poo i ncrea sed above 1. The curves have this characteristic shape 
even when plotted on a logarithmic scale for Pe/p^* 


In general, thrust coefficient decreased with increase of stagnation 
temperature T^ when the nozzles were overexpanded (pg/p^ < l), this ef- 
fect becoming more pronounced as either design pressure ratio or exit 
static pressure ratio is decreased. Only negligible effects of T^ on 

this potential-flow thrust ratio were found for underexpanded nozzles. 
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Above 4000° R stagnation temperature, the potential-flow thrust 
curves become more sensitive to P e /Poo as p^ is decreased. The most 

pronounced effect found is for a nozzle designed for a pressure ratio of 
5 when operated at an exit static pressure ratio of 0.24 and 6000° R 
stagnation temperature. Thrust ratio in this extreme case is reduced 
almost 20 percent as the stagnation pressure is reduced from 32 to 2 
atmospheres . 

Viscous flow . - Curves showing the viscous-thrust-loss parameter 
are also presented in figure 3. The ratio V F i is obtained after 

selecting and r*. Similar to the potential-flow thrust curves, 

these curves also show the greatest change with the variables considered 
when P e /P 00 is less than 1. The thrust-loss parameter always increased 

with decrease of exit static-pressure ratio. 

Optimum wall angle . - Figure 4 is presented to indicate how much 
thrust decreases with deviation of nozzle vail angle a from its value 
for maximum thrust; the carets mark the peaks of the curves. Thrust 
decreases more rapidly with decrease of a below an aerodynamic optimum 
than it does with increase of a above optimum. Only small differences 
in this rate of decrease with changes of stagnation pressure or temper- 
ature, wall temperature, design pressure ratio, or exit static-pressure 
ratio were found. 

The potential- flow thrust ratio (eq. (8)) is also shown in figure 
4 (dashed line) so that the amount of thrust loss due to exit flow angle 
and off-design expansion can be compared with that due to viscosity. At 
a Vq/Voo 1 j the two losses are the same* order of magnitude. At 
P e /Poo ~ 0.25 and * 2, both losses are large; but, at Pg/p^ 8 and 

dy * 20, the viscous losses are much smaller than that of the potential 
stream. 

Optimum wall angle in figure 5 is shown to decrease with increase 
of exit static-pressure ratio, design pressure ratio, stagnation pres- 
sure, and wall temperature and to increase with increase of stagnation 
temperature. Only sample plots are shown :.n figure 5, since optimum 
wall angle can be easily obtained by use of equation (l3) if figure 3 
is available for the desired conditions. 


Boundary-Layer Characteristics 

Momentum thickness . - Momentum thickness (eq. (3) ) was found to in- 
crease with increase of stagnation temperature and with decrease of wall 
temperature in the lover temperature range of 1000° to 3000° R (fig. 6). 
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In the upper temperature range , momentum thickness decreased with in- 
crease of stagnation temperature, wall temperature, or stagnation pres- 
sure. Values for the convergent part of the nozzle are shown only for 
the lowest stagnation pressure because of their crossing over the curves 
of the divergent part of the nozzle. 

Profile parameter . - Profile parameter (eq. (4)) is shown in figure 
7. The increase of profile parameter with increasing wall temperature 
and decreasing stagnation temperature is obvious from equation 4. No 
appreciable effect of stagnation pressure was found for 1000° = T^ 

= 4000° R. This effect of pressure appeared to be random and, since 
the scatter was less than 4 percent, only a mean line was shown in this 
temperature range. At temperatures much in excess of 4000° R, the pro- 
file parameter is consistently increased with increase of stagnation 
pressure. 


Heat Transfer 

The heat- transfer rate per unit area is shown in figure 8, and the 
heat-transfer rate through the whole nozzle is shown in figure 9. The 
curves are shown for the convergent as well as the divergent part of the 
nozzle. The maximum heat flux is shown to be at, or slightly downstream 
of, the nozzle throat. Both figures show the expected increase of heat 
transfer to the nozzle with increasing stagnation pressure or temperature 
and with decreasing wall temperature. Figures 8 and 9 are both for a 
constant wall temperature. Figure 10 presents the adiabatic wall tem- 
perature, that is, for the case of no heat transfer in either the x or 
y directions. The ratio of adiabatic wall temperature to stagnation 
temperature is seen to decrease with increase of stagnation temperature 
and with increase of stagnation pressure when Tq = 4000° R. 


Valid Range of Exit Static Pressure Ratio on Charts 

In practically all applications of an exit nozzle, a base region 
due to structural considerations will exist between the external air- 
stream and the exhaust jet from the nozzle. The pressure in this base 
region, p^, is the pressure to which the nozzle exhausts and can be, 
therefore, of prime importance in determining the internal flow of an 
overexpanded nozzle. The undisturbed external pressure p^ enters into 
the thrust equation by definition of jet thrust. When the nozzle is 
surrounded by still air, no difference exists between p^ and p^. 

When the nozzle is surrounded by a moving external airstream, however, 
as in flight, the base pressure is an involved function of the geometry 
of the vehicle encompassing the nozzle as well as a function of the noz- 
zle internal and external streams. This problem is discussed in refer- 
ence 3. 
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As operating pressure ratio p t /p ro is increased above the design 
pressure ratio (Pg/p^ = l) by decreasing ambient pressure p^, and hold- 
ing inlet stagnation pressure constant, the flow inside the now- 
underexpanded nozzle will remain the same, and no additional computational 
difficulties are encountered. The thrust ratio F/Fj varies, however, 
as the ideal thrust varies with operating pressure ratio or Pg/p^. 

Downstream of the nozzle exit, the flow will continue to expand until 
it reaches ambient pressure. 

As the operating pressure ratio is decreased below design pressure 
ratio by increasing p^ and holding p t constant, the flow inside the 

now-overexpanded nozzle will remain essentially the same as long as the 
shock that is then present stays at the exit station of the nozzle. 

Since the flow inside the nozzle then remains an accelerating flow, it 
tends to stabilize the boundary layer. The exit shock will not, there- 
fore, start to move into the nozzle until the nozzle operating pressure 
ratio is reduced low enough so that the pressure ratio P e / across 
the shock is large enough to separate the boundary layer. If the nozzle 
pressure ratio is decreased still further, the shock will move upstream 
into the nozzle and drastically alter the internal flow. At this point 
the charts of this report are no longer valid. When the shock is still 
at the exit station, any disturbance of th:.s shock on upstream conditions 
is considered to have a negligible effect on nozzle performances. The 
charts are therefore valid for 


Pe 



Theoretical and experimental evidence of reference 4 indicates that 
the pressure ratio required to separate an adiabatic turbulent boundary 
layer at a moderate temperature level must be at least as strong as the 
static pressure that occurs across a shock at a Mach number ratio of 
0.762. Reference 5 concludes that there if no effect of heat transfer 
on turbulent boundary-layer separation over the range of T w /T^ from 

0.75 to adiabatic. This result was obtained from schlieren photographs 
of the flow over a forward- facing step (coilar) mounted on an axisym- 
metric cooled body mounted in a supersonic wind tunnel. Additional data 
obtained at the Lewis Research Center from photographs of flow over a 
step mounted on a cooled cone indicated no effect of heat transfer even 
at a T w /T of 0.25 (fig. 11(a)). Figure ll(b) shows the corresponding 

separation pressure ratio for various Mach numbers ahead of the shock; 
this figure was computed from oblique shock charts and a constant Mach 
number ratio across the shocks of 0.762. 
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Only a value of p^/p^ is then required to estimate the pressure 
ratio at which a nozzle will first separate, since 

Pt Pb 

Poo Pe Poo 

(I t\ = = ~PtT 

\ P -A . p e 

7 design 

and 

Pb 

/Pe\ = Pc 

\Pco/ /Pb\ 

separation I ~ J 

\ p e/ 

separation 
SUMMARY OF RESULTS 

The thrust^ boundary-layer; and heat-transfer characteristics of a 
family of radial- flow nozzles were computed. Stagnation pressure was 
varied from 1 to 32 atmospheres; stagnation temperature from 1000° to 
6000° R; and wall temperature from 1000° to 3000° R. All calculations 
were made for air in chemical equilibrium. The results; where desired; 
were generalized independent of wall angle. 

The following results were obtained: 

1. There was no effect of stagnation pressure on profile parameter; 
adiabatic wall temperature; or potential- flow thrust ratio of the nozzles 
up to a stagnation temperature of 3200° R. Beyond this temperature; pro- 
file parameter decreased; adiabatic wall temperature increased; and the 
potential-flow thrust curves became more sensitive to off-design pressure 
ratios as stagnation pressure was decreased. 

2. Total thrust decreased more rapidly with decrease of wall angle 
below the value for maximum thrust than with increase of wall angle above 
that for maximum thrust. Optimum wall angle decreased with increase of 
exit static pressure ratio; design pressure ratio, stagnation pressure, 
and wall temperature; and it increased with increase of stagnation 
temperature . 

3. The maximum heat flux was near the nozzle throat. Both total 
heat transfer through the nozzle and heat transfer per unit area increased 
with increasing stagnation pressure or temperature and with decreasing 
wall temperature. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, October 15, 1958 
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APPENDIX A 

DERIVATION OF AVERAGE COSINE OF ANGLE BETWEEN POTENTIAL- 
FLOW VELOCITY VECTORS AND NOZZLE AXIS 


In radial flow, all flow properties are a function only of R, where 
R is the radius vector from the hypothetical point source to any spher- 
ical surface with the point source as center (fig- l(a)). Therefore, 
mass flow per unit area, as well as force per unit area, is constant at 
any given radius R if area is taken in one of the aforementioned spher- 
ical surfaces. Thrust in a direction parallel to the nozzle axis is 

/ A fk 

pU 2 cos a dA + / (p - p )cos a dA 

00 


0 


0 


where the symbols without subscripts are taken at the nozzle exit station. 
By definition, 


F = 


r T w tt 2 

2jt / 

"0 


cos a 


dr 


+ 2 * f 

0 


w (p - Poo)r 

■“55TST" djr l cos av a - 


i 


2k 


f 

0 


w 


2 r rw 

pU r dr + 2« / 

*0 


(p - P«>) r dr 


where 


, . 2itr 
dA = dr 

cos a 


Solving for cos av a and using the relation r = R sin a gives 


cos av a = 


/* r w /' a v 

j r dr / sin a cds a da 

_T) _ 

f Tv r dr / %a<w ’ i 

/ / sin x da 

J 0 cos a J Q 


1 + COS Oy. 


( 2 ) 
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APPENDIX B 


DETERMINATION OF PROFILE PARAMETER 


The compressible, adiabatic boundary-layer equations can be reduced 
to the incompressible form by use of a Stewart son- type transformation 
(ref. 6) . 



1 


Effects of compressibility can therefore be studied by use of these 
transformation equations (refs. 7 to 9). 


Reference 9 indicated that the incompressible profile parameter 
is constant at a flat-plate value in a nozzle with all favorable pressure 

gradient < 0^. In this reference, was obtained by operating on 

the measured compressible profile parameter H with a Stewartson trans- 
formation. It is assumed herein that the incompressible profile param- 

eter = -g- is constant at a value of 1.286, which corresponds to a 
y i 

l/7 -power profile, where 



When nonadiabatic flow such that there is heat transfer through the 
boundary layer is considered, an additional integral term arises in the 
transformation process. The purpose of this appendix is to discuss a 
basis for the evaluation of this term. The resulting equations for H 
are the same as (9) and (25) of reference 2 if is held constant, but 
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the derivation herein for this special case does not employ as in- 
volved or empirical a procedure. It is hoped that the following will 
add clarity to the assumptions involved. 

Displacement thickness for compressible flow is defined as 



by using the equation of state and assuming molecular weight and static 
pressure constant with respect to distance normal to the wall. Then, if 
specific heat is constant, 



where T is the stagnation temperature corresponding to the static tern- 
s 

perature T. Substituting equation (B5) in (B4) and making use of (Bl), 
as in reference 4, 



and 



where momentum thickness transforms as 


e = eM i 
HPp H 
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and the integral 


n - 


lldY remains to be evaluated. 


Incompressible turbulent boundary layers have velocity profiles that 


can be quite closely approximated by the relation 


(ref. 10). 


"p,i 


It can be shown (ref. ll) for a flat plate in incompressible flow 

U- n 

^ H i = constant) that, if is a so l u tion to the boundary- 

layer momentum equation, then ^ = (g~ ) is a solution of the 

energy equation when Prandtl number is 1. This follows from the fact 

T - T 

that the energy equation, when written in terms of © = =2 has the 

X - 

same form and boundary conditions as the momentum equation. Therefore, 


n / T w 


n + 1\T 


Also, 


Hi = 


§ /■©"[* - ffl'k 


n + 1 


[n + 1) ( 2n + I' 


= 2n + 1 


= 1.286 if n = y 


Substituting equations (B7) and (B8) into (B6) yields 


E -in MXT i + % + 


T t + 1.286T W ^ h t + 1.2861^ 


_i n U 
■* n + 1 VT 


l)H i + E 1 + 1 


1 


(4) 
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(b) Family of nozzles such that L2/L1 = constant. 


Figure 1 . - Description of nozzles. 
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Optimum wall angle, a w ^ 0 pt* de S 
Figure 2. - Plot of equation (13). 
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Potential-flow thrust, — — — — ^ Viscous thrust loss 












Potential-flow thrust, TTg — =r v*- Viscous thrust loss 









Potential-flow thrust, — — — — — E Viscous thrust loss 

























Potential-flow thrust, Viscous thrust loss 
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Potential-flow thrust, — - — — Viscous thrust loss 











Potential-flow thrust, — 7 ; tT Viscous thrust loss 
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Figure 3. - Continued. Thrust characteristics. 
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(a) Variation with exit static-pressure ratio and design pressure ratio. 
T w , 1000° R; T t , 1000° R; Pt /p std , 32; r*, 1. 


Figure 4. - Effect of nozzle wall angle on total thrust ratio. 
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Fi glare 4. - Concluded. Effect of nozzle wall angle on total thrust ratio. 


Optimum wall angle, S^opt’ de 6 
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Exit static-pressure ratio, p e /p 

(a) Variation with potential-stream stagnation pressure and exit static-pressure 
ratio. T w , 1000° R; T t , 1000° R; rj, 1. 


Figure 5. - Optimum nozzle wall angle. 




40 



L x> ‘ 3 "[Sub xx' bw tnrrarj (.dio 


Figure 5. - Concluded. Optimum nozzle wall angle. 
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(b) Effect of stagnation pressure. Stagnation temperature, 4000° R. 
Figure 6. - Continued. Momentum thickness. 
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(a) Stagnation temperature, 2000° R. 



Figure 7. - Boundary-layer profiles. 
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(c) Stagnation temperature, 6000° R. 

Figure 7. - Concluded. Boundary-layer profiles. 
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Figure 8. - Continued. Local heat-transfer rate. 
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(b) Stagnation temperalure, 4000° H. 


Figure 9. - Continued. Total heat 1 ransfer through nozzle. 


51 



50,000 


40,000 



10 12 14 16 

Area ratio, A/A* 

(c) Stagnation temperature, 6000° R. 

Figure 9. - Concluded. Total heat transfer through nozzle. 
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0 .4 .8 1.2 

Ratio of wall temperature to free- 
st ream stagnation temperature, T w /Trp 


(a) Separation Mach number ratio. 



(b) Separation pressure ratio. 

Figure 11. - Turbulent boundary- layer separation characteristics. 
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